Energy landscapes in random systems, driven interfaces and wetting 
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We discuss the zero-temperature susceptibility of elastic manifolds with quenched randomness. It 
diverges with system size due to low-lying local minima. The distribution of energy gaps is deduced 
to be constant in the limit of vanishing gaps by comparing numerics with a probabilistic argument. 
The typical manifold response arises from a level-crossing phenomenon and implies that wetting in 
random systems begins with a discrete transition. The associated "jump field" scales as (h) ~ L~^^^ 
and L~^'^ for (H-l) and (2-|-l) dimensional manifolds with random bond disorder. 
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The physics of systems with quenched disorder is 
related to the energy landscape. The free energy is 
at low temperatures governed by zero temperature ef- 
fects, which in turn are ruled by the scahng of the 
disorder-dependent contribution. Random magnets, as 
spin glasses and random field systems, flux line lattices 
in superconductors, and granular materials are exam- 
ples of physical systems in which frustration and disorder 
play an important role. Disorder may dominate also in 
non-equilibrium conditions, like driven systems (domain 
walls in magnets, flux lines in superconducting materi- 
als). In that case temperature-driven dynamics (creep, 
aging) and the external drive change the system from one 
metastable state to another 0J^. 

A lot of information about energy landscapes is con- 
tained in how the number of local energy minima and the 
typical scale of their energy differences scale with system 
size, L This can be interpreted in a geometric fashion 
in that one compares the energy difference of two states 
with their overlap in terms of the spin configuration (as 
for magnets). In spin glasses an intense debate still goes 
on: whether in the thermodynamic limit the thermody- 
namic state is trivial ( "droplet" picture B) or not (as in 
the "repHca symmetry breaking" picture^^). 

Consider now the problem of the energetics of D di- 
mensional elastic manifolds in random media (|-|^, of 
which the best-known case is a directed polymer (DP) 
in a random medium with D — 1, often called a 'baby 
spin-glass' |l^. For these systems the interface energy is 
proportional to the area, and the sample-to-sample en- 
ergy fluctuations scale with the exponent 6, {0 = 1/3 for 
a DP m d — D + 1 = 2 embedding dimensions). The ge- 
ometry is often self-affine, characterized by a roughness 
exponent (2/3, when d = 2). In the simplest energy 
landscape the valleys and excitations are separated by 
energy gaps proportional to where / is the length scale 
of the perturbation. 

Here the susceptibility of elastic manifolds is studied 
in the presence of weak fields numerically and by scal- 
ing arguments. By investigating each sample separately, 
we explore the changes in the energy landscape with ap- 
plied fields. These lead to discrete 'jumps' in the phys- 
ical configuration. As a consequence scaling arguments 



of wetting in random systems do not work in the limit 
of weak fields if the original interface-to-wall distance is 
much larger than the interface roughness [|ll| . With pre- 
conditioned systems we obtain the detailed probability 
distribution of the energy differences (gaps) between lo- 
cal minima and the global one. We find that the average 
interface behavior can be explained with scaling argu- 
ments, but the susceptibility can not, and it is directly 
related to the exact properties of the gap distribution. 
Thus the detailed statistics of the landscape is important. 
This contradicts considerations for random systems that 
assume well-defined thermodynamic functions [ p^ and 
scaling arguments with a single parameter (L^). These 
findings agree with claims that the susceptibility of a DP 
to thermal perturbations or applied fields, is anomalous 
p3|-p^. The reason is that the response to a very weak 
field, say applied locally at the end-point of a DP, is gov- 
erned by rare samples. The disorder-averaged response 
differs from the typical one because the ground state can 
be almost degenerate with a local minimum. Likewise, 
numerical studies of c? = (1-1-1) DP susceptibility reveal 
aging phenomena reminiscent of real spin-glasses |p^ , |l7| . 

The continuum Hamiltonian for a D dimensional elas- 
tic manifold (x is an internal coordinate and a z (scalar) 
displacement) 

H = y d^X [r{V2(x)}2 + Vrix, Z) + hiz)] , (1) 

with an elastic energy (F is the interface stiffness), and 
Vr a random pinning energy (we use a random bond cor- 
relator, {Vr{-X,z)Vrix',z')) = 2VS{x - K')d{z - z')) . h{z) 

couples the interface to an external perturbation, e.g. it 
describes a constant magnetic field H in Ising magnets 
with antiperiodic boundary conditions. 

The Hamiltonian (0) describes also complete wetting in 
a random system, where h{z) equals to the chemical po- 
tential difference of the wetting layer and the bulk phase 
]ll| , p^ [2C|] . For h non-negligible the wetting- inducing ex- 
ternal potential competes with the tendency of the in- 
terface to win pinning energy. Assuming that these bal- 
ance, the average interface-wall separation (z) becomes 
(z) ^ h~'^,ip — ^^T^-j;: where "ip is the depinning exponent. 
T measures the scaling of the elastic and pinning energy 
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and is given by t = 2(1 — and n is the scaling expo- 
nent of the external field h{z) ^ (here we use k = 1). 
For random bond systems T~lmd=l + l dimensions, 
and r ~ 2.9 in d = 2 + 1 using the known bulk roughness 
exponent values 2/3 and 0.41 in d = 2 and 3, respectively 
P,pil. In d = 2 numerical simulations in random Ising 
systems indicate, in agreement, -0 ~ 0.5 [ p^|JlS| ]. 

A network flow algorithm, invented by Goldberg and 
Tar j an |^^, is used here for the numerical procedure. It 
solves the minimum cut - maximum network flow prob- 
lem, and produces in polynomial time the exact ground 
state energy and interface configuration given a sample 
{L X Lz 01 L X L X Lz) with fixed quenched disorder. Lz 
is the z-directional system size. The algorithm is con- 
venient when one makes systematic perturbations to the 
original problem [h — 0) [g3[|4j. Figure |l| illustrates the 
sample-to-sample behavior, as the external field h{z) is 
switched on slowly (see Eq. (||)). At /i = the interface 
is in the ground state. It has a mean wall distance zq and 
a width w ^ in a system of transverse size Lz ■ As the 
field is increased the interfaces move intermittently with 
jumps to positions (zi, Z2, ■ ■ ■ , z„, . . .) This corre- 

sponds to a first-order transition. Instead of finite-size 
excitations the first change in the interface configuration 
is a macroscopic jump with zero overlap between the old 
and new states. The first transition point defines a jump 
field hi. It assumes the role of a latent heat, and corre- 
sponds to the landscape-dependent energy to move the 
interface. 

The two possible mechanisms are compared in the inset 
of Fig. |l]. Either the interface adjusts itself gradually by 
forming 'bubbles' or local excitations, or it jumps com- 
pletely (compare with the main figure). The scenarios 
are linked to the structure of the energy landscape. If 
the first excitation is localized and has the transverse 
spatial extension A (Z ~ l^^/^) [|^, the energy cost 
scales with A"/*" and the energy win in the field scales 
with /ii A^+'^''^^^/''. Assuming that a — 6 the jump field 
hi ^ /S.^ = /S.'^/i-^-id--^)/^. The exponent is negative, 
and thus small excitations are the more expensive ones 
p6| . Numerically, the fraction of jumps leading to a non- 
zero overlap with the ground state decreases towards zero 
slowly with L. Also, the scaling function of the interface 
jump lengths approaches a constant shape. The mean 
jump length (Azi = — zi, zi < zq) scales extensively, 
Azi ~ Lzj not with e.g. L'^ . 

So for small fields h and <C Lz the sample-to-sample 
fluctuations lead to a discrete (wetting) transition. The 
average behavior with {z{h)) and typical interface behav- 
ior with z{h) do not coincide, since the asymptotic — > 
limit is dominated by the near-degeneracy of the ground 
state. In the limit L'-' <^ Lz there are many independent 
'valleys' in the energy landscape for directed surfaces. 
Each of these has an energy En corresponding to a local 
minimum and their energy difference to the ground state 
(with Eq) is expected to scale as with two independent 
sets of disorder. That is En ~ Eq ^ L^ . This energy 



difference equated with the jump energy /iiL^Azi leads 
(with the choice Lz = L) to the scaling 

/ii ^ L^-'^ = i-". (2) 

The jump field exponents are a = 5/3 and a 2± 2.18 
in d = 2 and d = 3 random bond systems, respec- 
tively. In d = 3 random field interfaces have a = 5/3 
(C = 2/3 and 6* = 2C £> - 2 §||). It is assumed 
that Azi ~ i, since the valley energies are indepen- 
dent, except for the bias caused by the field h. Figure 
H compares the exponent values to numerical data with 
only the non-overlapping jumps being considered (with- 
out this pruning the same exponent is obtained asymp- 
totically). For D = 1 a becomes 1.62 ±0.04, close to the 
scaling estimate of 5/3. The inset shows the disorder- 
averaged jump distance (Azi) vs. L and shows that the 
interface response geometry scales linearly with L (as 
discussed above). For D = 2 random bond manifolds 
we obtain a ~ 2.2, in reasonable agreement again. In 
the limit {zn(h)) ~ z„(/i) 'i^ w ^ L'^ (after n jumps of 
sizes Az„ = Zn-i — Zn) the mean- field wetting theory ap- 
plies, and indeed we obtain for the depinning exponent 
for d = 2 ^/j ~ 1/2, and for d = (2 1) i/" ~ 0.26, in rough 
accordance with the Lipowsky-Fisher [ pT| prediction. In 
d = (2+1) there are deviations including a dewetting 
transition for weak disorder | |2l[ | and the exponent con- 
verges very slowly (zq ~ w at L ~ 10^ if Lz — 50 

. . . . 

If the initial interface position is random, the jump 
statistics are an average over the initial number of avail- 
able valleys (recall that the field breaks the up-and-down- 
symmetry, see Fig. |l|). Thus we also consider the limit 
in which the initial position is set to be inside a fixed- 
size window, Zq/ Lz — const. We expect that the number 
of local valleys in the landscape, accessible with /i > 0, 
has a well-defined average (in the grand-canonical sense), 
and that the relevant scaling parameter is Lz/L'^ . Fig- 
ure ^ shows the scaling function of the probability distri- 
bution P{hi) obtained with this initial condition. We 
find the form P{hi/{hi)) = A{L) f {hi / (hi)) where A 
depends on the energy gap scale L^ and / is a scaling 
function with the limiting behaviors f{x) ^ l,a; — > 
and f{x) ~ exp {—ax^),x > l,/3 ~ 1.3. The distribution 
is constant for small fields and has an almost exponen- 
tial cut-off. The scaling properties imply in particular 
that the disorder-averaged susceptibility diverges. The 
change in magnetization is given by the number of in- 
terfaces that have moved times the mean distance (Azi). 
Thus the divergence is not Xtot ~ L^ [|l^ . Figure ^ shows 
the average jump field in the fixed height ensemble with 
varying Lz and constant L. We have fitted the data with 
(hi) ~ L~'^ , and the best fit is obtained by the scaling 
exponent 7 ~ 4/3. 

Consider now the energy landscape for small h. It has 
k = 1, . . . ,Nz associated minima {Nz ^ Lz/L'^) with the 
energies Ek picked out of an associated energy gap prob- 
abihty distribution P{AEk), where AEk — Ek — Eq and 
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Eq is tlie ground state energy. When /i > 0, all the local 
minima attain an energy of Ej- + hAzk with respect to 
the reference state with zq and Eq. Now we make the as- 
sumption, analogous to the Random Energy Model , 
that all the gap energies AEk are independent random 
variables. We can now simply compute the probability 
for the original ground state being stable for any h (i.e. 
no jump has taken place) by the joint probability Pq that 
all the Ek + hAzkS are still higher than the original one 
with the given h. —dPg/dh gives then the probability 
that this level crossing occurs at exactly h. By comput- 
ing 

dl^^_ - J"' p'"^ Pi.)dxdk P{kh/N,) 
dh ' ' " A 1^1^''/''' Fix) dx 

(3) 

one can show that the only P that reproduces the numer- 
ical P{hi) is a constant one, whereas all other functional 
forms of P fail, see Fig. ^. This P is in fact exactly 
the marginal one needed for the susceptibility per spin 
X ~ limh^Q (dz/dh) to diverge in the thermodynamic 
limit. In particular for a distribution P{hi) that van- 
ishes in the zero field limit the susceptibility would stay 
finite. Using the obtained form for the probability distri- 
bution gives X ~ i jjY "^here 7 ~ 1 — ^ relates to the 
density of valleys. This slightly disagrees with the above 
result (7 c± 4/3) since with L = const x ^ ^2, — 1- 
In the isotropic limit L cx the extensive susceptibility 
simply reads Xtot = L\ ~ L''+i+^-^ ~ 1"^^^+^. To con- 
clude, X (or Xtot) is determined by the exact low-energy 
properties of P, or by the rare events in the low AE tail. 

To summarize we have studied the coupling between 
the energy landscape structure and the response of in- 
terfaces, related for instance to complete wetting. A dis- 
order averaging that reflects correctly the level-crossing 
character of the problem reveals that the wetting starts 
with a discrete transition. Thus the randomness of the 
energy landscape drives a second-order transition to a 
first-order one. The 'jump' is associated with an effective 
specific heat, which can be understood in terms of scaling 
arguments. The susceptibility is governed by the infre- 
quent cases with low-lying local minima, which allows us 
to derive a constant energy gap probability distribution. 
The results should be relevant for other problems like flux 
line lattices in superconducting materials with quenched 
randomness It will also be of interest to see if the 
energetics and the geometrical character of the response 
can be coupled with arguments concerning the energy 
barriers in each specific configuration |p9| . This would 
allow to understand the dynamics in the creep regime, 
when the interface moves between metastable states. 
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FIG. 2. Finite size scaling of the average first jump field 
(hi) for one dimensional DP's. The line is the least squares 
fit to data. The scaling argument gives a = 5/3. The inset 
shows the average jump distance (A^i) at the corresponding 
field hi with a linear fit to data. { ) is the disorder-average over 
A'^ — 1000 realizations for the system sizes L x Lz = = 50^ 
and 100^ iV = 500 for = 200^ ~ 400^ and N = 200 for 
= 600^ — 1000^ . The disorder is of random bond type. 



N 

c 

o 



o 
a, 



13 



1.0 
0.8 
0.6 
0.4 
0.2 



0.0 



1 , z„(x) 



100 - 



.S." 50 



1 , z,(x) 
"2", zjx) I 



50 100 150 200 
X 



lO"' 10 10 10 10" 

applied field (h) 



10" 



FIG. 1. Overview of two realizations of changes in mean 
heights z of interfaces normalized by their original (global 
minimum) positions zo vs. applied field h for (1 -f 1) di- 
mensional systems. Note the large jumps in both cases. 

= 200^. Jij,z £ [0, 1] uniform distribution and Jij^x ~ 0.5 
(random bond disorder). The expected scenarios (bubble for- 
mation, jump to the lower edge of the system) before and 
after the first moves from global minima zo{x) to zi{x) are 
shown in the inset. 
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FIG. 3. The scaling function of the probability distribution 
P{hi/ (hi)) X {hi} for the first jump field values hi normalized 
by their disorder-average (hi) in a (10-base) semilog-scale for 
the system sizes L x Lz = L'^ = 100^ and 200^ The in- 
set shows the tails in the natural-log-scale. The initial global 
minimum position zq/Lz — const for all L. The number of 
realizations N = 10^ for both system sizes. The line is the 
analytic result from Eq. (^) with a uniform distribution P{x) 
and Nz =20. 
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FIG. 4. The disorder-average of the first jump field (hi) as 
a function of transverse system size Lz for the system sizes 

L =100, 150, 200, 250 and 300, each with zo/L^ = const. 
The number of realizations ranges from A'^ = 500 for L = 300, 
Lz = 500 to AT = 2600 for L = 200, = 600. The line 
L~'^,'y = 4/3 is a guide to the eye. 
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